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We discuss the rough surface effects on a two-dimensional chiral k x + ik y superconductor. 
The atomic scale roughness at the surface is considered using the random 5 matrix model. 
The roughness effects on the self-consistent order parameter, the surface mass current and the 
surface density of states are studied using the quasi-classical theory. We find that the surface 
mass current is suppressed by the surface roughness. The surface density of states shows a 
quite similar behavior to that of superfluid 3 He B phase. When the surface is specular, the 
surface Andreev bound states form a band which fills the bulk energy gap Abuik- When the 
surface becomes diffusive, there occurs a sharp upper edge of the surface bound states band and 
there opens a sub-gap between the edge and the bulk energy gap. We show that this sub-gap 
is induced by the repulsion between the surface bound states and the propagating Bogoliubov 
quasi-particles through the second order process of roughness. 
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The 2-dimensional chiral k x + ik y state is known as 
a model system for S^RuC^ superconducting state. 1 ) 
When there is a surface along the y axis, there occur gap- 
less surface Andreev bound states as in other p-wave 2, 3 ) 
and (i-wave 4 ) pairing systems, because the k x component 
of the order parameter changes its sign under the surface 
reflection. The surface bound states in chiral system are 
known to carry spontaneous mass flow along the sur- 
face. 5 ' 6 ) Recently, surface bound states are recognized as 
edge states which reflect the topological nature of the 
bulk pairing state. A lot of attention has been paid to 
the surface Andreev bound states from this aspect. 6 ~ 15 ) 

In this paper, we consider the effects of atomic scale 
surface roughness on the chiral k x +ik y state. We use the 
quasi-classical theory 16, 17 ' developed for the study of p- 
wave Fermi superfluids. We calculate the self-consistent 
order parameter, the surface mass current and the sur- 
face density of states. The surface density of states shows 
quite a similar behavior to that in the 3-dimensional BW 
state. 17 ) Existence of the order parameter component 
parallel to the surface disperses the surface bound state 
energy. The bound states, therefore, form a band below 
the bulk energy gap Abuik- When the surface is specu- 
lar, the band completely fills the bulk energy gap. When 
the surface becomes diffusive, however, there occurs a 
sharp upper edge of the band and there opens a sub-gap 
between the edge and the bulk energy gap. The band 
edge energy A* increases as the roughness is reduced. 
Similar sub-gap has been known in the BW state. 17 ~ 19 ) 
The sub-gap was first reported by Zhang 18 ) who treated 
the surface roughness using the thin dirty layer model. 
He suggested that the sub-gap is due to the suppres- 
sion of the parallel component of the order parameter 
by the roughness. However, Nagato et al. 17 ) found that 
this sub-gap also occurs when the order parameter is as- 
sumed to be spatially constant. Although the existence 
of the sub-gap played a decisive role in the interpretation 
of the transverse acoustic impedance of the B phase of 



superfluid 3 He, 20 ~ 23 ) the origin of the sub-gap has been 
a puzzle for a long time. 

We consider a two-dimensional k x +ik y superconductor 
which fills the x > domain. The surface extends along 
the y axis. Since the order parameter is suppressed near 
the surface, the order parameter will take a form 

A(k,x) = A±(x)k x +iAn(x)k y = A±(x) cos 0+iAii (x) sin <j>, 

(1) 

where (j> is the angle between the Fermi momentum and 
the x axis. 

We investigate the effects by surface roughness of 
atomic scale using random S-matrix model. 16 ' 17 ' 23 ) The 
surface is characterized by an S-matrix for the quasi- 
particles at the Fermi level in the normal state 
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where k y (q y ) is the y component of the incident (scat- 
tered) Fermi momentum and rj is a Hermite matrix that 
specifics the surface roughness. We assume that 77 is a 
random Hermite matrix which obeys rjk y q y = and 



with W a pa- 



rameter that specifies the roughness of the surface. One 
can show that W = 1 corresponds to the diffusive sur- 
face boundary condition and W = corresponds to the 
specular surface boundary condition. 16,17 ) 

Taking into account the surface roughness within the 
self-consistent Born approximation, we obtain the quasi- 
classical Green's function at the surface. 17,23 ) 

G ± (k y ,Q)=Gs + (G s ±t) „J - (GsTi), (3) 
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Here G+(k y ,0) and G-(k y ,0) are the quasi-classical 
Green's function for the Fermi momentum (k x > 0, k y ) 
and (k x < 0,k y ), respectively. G$ is the quasi-classical 
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Fig. 1 
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(color online) Self-consistent order parameter of k x + ik y 
state is plotted against the distance from the surface scaled by 
the coherence length £ = vp / irT c for some typical values of the 
roughness parameter W . 
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Fig. 2. (color online) Temperature dependence of the total cur- 
rent J y along the y axis. Solid curve is for the specular sur- 
face (W = 0) and the dotted curve is for the diffusive surface 
(W = l). 



Green's function at x = for the specular surface and £ 
is the surface self energy which is induced by the rough- 
ness. The angle average in the two-dimensional system 
is 

f.TT/2 

d(j) cos <j> ■■■ . (5) 
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The quasi-classical Green's function G±(k y ,x) at finite 
x is calculated by evolution operator technique. 16,17 ) 

Using the Matsubara Green's function, we can cal- 
culate the self-consistent order parameter and the edge 
mass current. In Fig. 1, we show the self-consistent order 
parameter at T = 0.2T C . Since the bulk energy gap is 
isotropic in the two-dimensional k x + ik y state, the or- 
der parameter shows a quite similar profile to that of the 
three-dimensional BW state. The perpendicular compo- 
nent A±(x) is suppressed near the surface. In case of the 
specular surface, 5 ) the parallel component Ay (a;) is en- 
hanced near the surface such that compensates the loss 
of the condensation energy caused by the suppression of 
the perpendicular component A±(x). In case of the dif- 
fusive surface, Ay (x) is also suppressed by the incoherent 
phase mixing during the reflection processes. 24 ' 

Once the order parameter is determined, the surface 
mass current along the y axis can be calculated from the 
diagonal element of the quasi-classical Green's function. 
In Fig. 2, we show the total surface mass current J y , 
current density integrated over i, as a function of tem- 
perature. In case of the specular surface (W = 0), the 
total current tends to J y = —nh/4 as T — > OK. 6 ' Here, n 
is the total number density. When the surface is diffusive 
(W = 1), the total current is definitely suppressed. The 
suppression of the mass current by surface roughness was 
discussed by Ashby and Kallin 24 ) using GL theory. 

The surface density of states can be calculated from 
the quasi-classical Green's function with real frequency 
e. We show the angle resolved density of states in Fig. 3 
for the diffusive surface (W = 1). In case of the specular 
surface, the surface density of states shows a delta func- 
tion peak that corresponds to the surface Andreev bound 
state. The peak position is roughly equal to Am (0) sin^>. 



When integrated over the angle, therefore, the bulk en- 
ergy gap below Abuik is filled by the bound states. In case 
of the diffusive surface, the bound state peak is broad- 
ened and is shifted towards the lower energy. Moreover, 
there appears a sharp upper energy edge A* common to 
all the incident angles, which leads to a sub-gap between 
A* and A bu ik- 
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Fig. 3. (color online) Angle resolved surface density of states in 
case of the diffusive surface (W = 1). Incident angles are = 

°> i". M*" and IK 



To examine the origin of the sub-gap, we consider 
Green's function at the surface given by Eqs. (3) and 
(4). To discuss the density of states below the bulk en- 
ergy gap Abuik, we consider real frequency |e| < Abuik- 
Let us first consider Gs which is given by 17 ' 23 - 1 

^(1 + L> 2 ) -2D 

-i{l + D 2 



G, 
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(6) 



where D = D(0, e, k y ) is a solution at x = of the Ricatti 
equation 

d „ 
v F cos (f>— D = -2ieD + A(fc, x)D 2 - A*(k, x) (7) 
ox 
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with the boundary condition at the bulk infinity 

iA*(fc,oo) 



D(oo, e, k y ) 



J(a-!p) 



(8) 



e+V^-IAbuikl 2 

where we have defined a = sin _1 (e/Ab u ik)- For the en- 
ergy |e| < Abuik, it can be shown from Eq. (7) that 
\D(x, e, k v )\ is always unity, therefore we may write 



D(x, e, k y ) 



(9) 



with 6 the real function. Solving Eq. (7) for a given k y , 
we find an e that satisfies 



D(0,e,fc y ) = 1. 



(10) 



This energy is the surface bound state energy for the 
specular surface because the Green's function Gs has a 
pole at that energy. It is worth noting that D is related 
to the Nambu amplitude (u(x),v(x)) of the state with 
energy e. The ratio of the hole component v to the par- 
ticle component u is given by v/u = {—i)D. It follows 
that the hole-particle ratio v/u at the surface is equal to 
—i for all the surface bound states. 

When we assume that the order parameters are con- 
stant, i.e., A±(x) = A||(x) = Abuik, D is also a constant 
given by Eq. (8). The bound states have a linear dis- 
persion relation e = Abuik sine/) = Abuik^y and can be 
regarded as Majorana-Weyl Fermions. 6 ) 

Now we consider the surface self energy. From Eqs. (4) 
and (6), we may parametrize the self energy in a form 



(11) 



is 3 Si 
si -is 3/ 

It is convenient to introduce projection operators P± = 
^(1 ± P2) with p2 a Pauli matrix in particle-hole space. 
Then we can write 

£ = ip 3 (P+(s 3 + Sl ) + P_(s 3 - si)) , (12) 



G s = i P3 P 



1-D 
1 + D 



1 + D 
1-D 



(13) 



It is obvious from Eq. (13) that P+ projects out the sur- 
face bound states. From Eqs. (4), (12) and (13), we find 
that 



s 3 + si = 2W 



s 3 - Sl = 2W 



1 - D 



(1 + (s 3 - Si)) + (1 - (S3 - 81)) D 



(14) 



1 + D 



(1 + (S3 + Sl))-(1-(S 3 + Sl))£> 



(15) 

The density of states is given by the imaginary part of 
the diagonal element of the quasi-classical Green's func- 
tion given by Eqs. (3), (12) and (13). 17 > 

^11 



G L + L =Gi 



i_ / (1 + (s 3 + si)) + (1 — (s 3 + si)) D 
'2 ^(l + (s 3 + si))-(l-(s3 + si))L> 

(1 + (s 3 - Si)) - (1 - (* 3 - Si)) D 



(l + (s 3 -si)) + (l-(s 3 -s 1 ))£>y (16) 
When the energy is in the range |e| < Abuik, D is given 



from Eq. (9) by a form e te {o,e,k v )^ j t f n ows that if both 
s 3 and si are pure imaginary, the diagonal element of 
the Green's function is real, namely there is no density 
of states. At the sub-gap energies, therefore, S3, si are ex- 
pected to take pure imaginary values. At first sight, both 
Eqs. (14) and (15) have pure imaginary solutions. Both 
the equations are invariant under the complex conjugate 
transformation because D = e 10 . The real part emerges 
when there appears a pole along the angle integral in 
Eqs. (14) and (15). 

From now on, for simplicity, we consider a case where 
the order parameters are constant and the roughness 
parameter W is small. When the order parameter is 
constant, D — e ! ( Q_< ^; therefore, the bound state en- 
ergy for the specular surface is given from D = 1 by 
a = 4> (e = Abuik sin </>) . Within the lowest order correc- 
tion with respect to W, s 3 + si remains pure imaginary 
because the bound states are projected out in Eq. (14). 



S3 
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On the other hand, S3 — Si acquires real part because 
Eq. (15) has a pole of the surface bound state when S3+S1 
is neglected. Let us consider the next order correction by 
W to the possible pole of Eq. (15) 

D - i±fa±ill . (19) 

When W is small, the possible pole will occur near D ~ 
1 + i(a — 4>). Expanding both sides of Eq. (19) in terms 
of small quantities, we obtain 

= (-2i)(s 3 + «i) = /(a). 



a 



(20) 
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Fig. 4. (color online) a 
W = 0.1 are plotted. 



- tt/2, tt/3, tt/2) and /(a) with 



We plot the both hand sides of Eq. (20) in Fig. 4 as 
functions of a — sin _1 (e/ Abuik)- Since f(a) is a decreas- 
ing odd function of a, there is a solution of Eq. (20) for 
any —tt/2 < <fi < tt/2. But, when we define a* at which 
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ir/2 and f(a) crosses; 

♦-7T/2 = /(<**), 



(21) 



we find that there is no solution of Eq. (20) for a in the 
range a* < \a\ < ir/2. It means that Eq.(15) has no pole 
and S3 — S1 remains pure imaginary in that energy range. 
As a result, for all the incident angles there occurs a com- 
mon sub-gap in the energy range A* = Abuik sm a * < 
|e| < Abuiio as seen m Fig. 3. Solving Eq. (21) with re- 
spect to a*, we obtain A* = Abuik sin a* as a function 
of W. The result is plotted in Fig. 5 together with the 
self-consistent solution of Eq. (4). 
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Fig. 5. (color online) A* vsW. Order parameters are assumed to 
be spatially constant. Dashed curve is a result of self-consistent 
solution of Eq. (4). Solid curve is a result of Eq. (21). 



The origin of the sub-gap is interpreted in a following 
way. Since Eq. (19) is also a possible pole of the diago- 
nal element G^. 1 of the Green's function (see Eq. (16)), 
Eq. (20) is interpreted to be an equation to determine 
the energy of the bound state with finite W, although 
we have used it to find out the energy range without 
solution. Equation (20) with (17) has a similar form to 
the denominator of the usual Green's function for the 
impurity problem within the Born approximation, there- 
fore it corresponds to the Brillouin-Wigner perturbation 
formula 



c(0) 



E 



IK, 



(22) 



where e„ corresponds to Abuik sin and |Kim| 2 to W. 
Since the bound states are projected out in Eqs. (17) and 
(20), the intermediate states are the propagating Bogoli- 
ubov quasi-particle states with energy |e| > Abuik- The 
right hand side of Eq. (22) becomes a decreasing odd 
function of e and reproduces the a dependence of f(a). 
The sub-gap comes out, thus, as a result of the repul- 
sion between the bound state and the propagating states 
through the second order process. This scenario does not 
change in case of the self-consistent order parameter, al- 



though Eq. (20) should be calculated numerically. The 
sub-gap in supcrfluid 3 Hc-B can be explained in a simi- 
lar manner. 

It is of interest if the sub-gap which has been observed 
in the B phase of superfluid 3 He can be also observed in 
Sr2Ru04, for example by tunneling experiment. 25, 26 - ) For 
comparison with experiment, the effects by finite trans- 
mittance of the rough interface should be examined. 27 ^ 
Such a study shall be reported elsewhere. 
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